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c3 ■ Abstract 

Penguin diagrams for decays such as 6 — > (s, (i)7 involve virtual loops of u or 
other light quarks. Logarithms of the virtual quark mass could, in principle, 
influence the phenomenological analysis of the decay. It is thus important 
to study these logarithms to all orders in QCD perturbation theory. In this 
paper we show that, at arbitrary order, the matrix elements of operators in 
the effective hamiltonian contributing to 6 — *> S7 are finite for the limit of 
m^ — > in penguin loops. 
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1 Introduction 

In the spectator- quark model, the inclusive weak decays of the 5-nieson 
can be pictured as QCD corrected 6-quark decays. Since in perturbative 
QCD, virtual light quarks appear in loop diagrams, their presence might 
suggest an enhancement of the rate through terms involving powers of In mioop 
(we denote by mioop the mass of any hght quark circulating in the penguin 
loops). The possible effects of such terms "long-distance" contributions were 
explored, for example, in Ref. |]I|. On the other hand, the finiteness for 
fnioop —>■ has been invoked in other papers, such as @, |3, §]. In this paper, 
we address this problem in the context of effective field theory. Our analysis 
is perturbative and therefore relates to current, rather than constituent quark 
masses |^. 

We refer to a particular flavor changing neutral current process, 6 ^ 57 
decay, but our results can easily be extended to other interesting FCNC 
processes, such as 6 — >• dj. We argue that, to any order in perturbative 
QCD, the limit mioop ^ in gauge invariant penguin amplitudes is finite, 
and the presence of virtual light quarks in these internal loops does not 
result in logarithmic divergences as mioop -^0. To be specific we argue that 
powers of In m^ioop are always accompanied by positive powers of m,ioop [§] • For 
example, in 6 — >■ (s, d)^ decays there is a contribution from penguin loops 
with a virtual u quark (jnioop = m^u), but the amplitude does not diverge as 
m,u — ^ 0. For definiteness, we discuss the exclusive b ^ s'j final state. 

Of course, penguin loops are not the only source of dependence on light 
quark masses. In the partonic 6 — > S7 amplitude, subdiagrams involving glu- 
ons (photons) attached to the outgoing s-quark as well as soft gluons attach- 
ing the 6-quark with the s-quark jet, are collinear divergent in perturbation 
theory in general. In practical calculations, these regions give logarithms of 
m,s that are not suppressed by powers of m^. In addition, the amplitude for 
any exclusive final state contains purely infrared divergences associated with 
the masslessness of the gluon. These kinds of divergence have already been 
treated in the hterature (see for instance @, 0, |^). As part of our analysis 
on penguin loops, we shall show that no additional logarithms arise from 
light quark loops collinear to the outgoing photon in the particular case of 
the two-body b ^ s'j decay amplitude. We derive our results by building 
the effective hamiltonian step by step, showing at each stage which kinds of 
dependence on light masses are to be expected, and by applying analyticity 
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arguments and IR power counting techniques developed in perturbative QCD 



2 The Effective Field Theory 

The effective field theory for i?-decays describes the physics at the scale 
H c::i mij <^ m^Y, after the top quark and the W boson have been integrated 
out. The effective hamiltonian at the lowest order in Oem is defined by a 
sum of local operators whose matrix elements between initial and final states 
reproduce the amplitude at low energies 
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where Vckm denotes the appropriate factor (typically quadratic in CKM 
matrix elements). The utility of the effective hamiltonian formalism is that 
it separates short distance contributions, described by the coefficients, which 
can be calculated perturbatively, from long distance contributions, incorpo- 
rated in the matrix elements of the local operators. 

We can summarize the steps to build the effective hamiltonian for h — >■ 
57 as: (i) calculating the coefficients by matching the full theory onto the 
effective theory at high scale; (ii) evolving the coefficients down to the scale 
0{mb) by the renormalization group (RG) equations; (iii) evaluating the 



matrix elements of the operators. In the same spirit of Ref. |T2], we will 
analyze the role of the light masses during these three steps. 

2.1 Matching 

The first step consists in matching the effective theory onto the full theory. 
To match means to extract the coefficients Cj by comparing the amplitude 
in the full theory and in the effective theory at the same order in a^- At the 
matching scale, all IR behavior cancels and logarithms of light masses are 
then eliminated in the coefficients. 

As an example, we shall discuss weak penguin diagrams in 6 — > 37, like the 
one shown in Fig. 1. The effective flavor changing gauge invariant couplings of 
the photon are of the type F^yd^'SL'~i^bL and 'SLF^^a^^h^. In a renormalizable 
theory like the standard model, the corresponding amplitudes must be UV 



finite, because the above operators have dimension higher than 4 and cannot 
arise as counterterms. The aniphtude of the full theory, to which the effective 
theory is matched, can be calculated by expanding in the ratio g/miy, where q 
is the photon momentum. This expansion can introduce infrared divergences. 
Since the external photon is taken on shell, the mass of the quark in the loop 
of Fig. 1 acts as an IR regulator. 

The expansion in photon momentum results in terms of the form Fi {q'^'j'^ 
— g^7 ■ q) and F2 ia^^q^. The form factor Fi includes an IR-sensitive term 
"^{xi — 1) InXj, where Xi = mf/m^^, with rrii the mass of quark i, circulating 
in the loop, while in F2 Inxi appears only multiplied by powers of Xi (see, 
for instance. Appendix B of Ref [|lll).[| 

In the effective hamiltonian, the local operators contain only light quark 
fields; the heavy quark fields have been integrated out. In the corresponding 
diagrams for 6 — i> 57, there is an analogous term of the type |ln(r7ij^//i^), 
where /z is introduced to fix the scale where the operators are renormalized, 
through the introduction of a counterterm proportional to F^^d'^'SL'^^hL. The 
coefficients of the effective hamiltonian are found precisely from the difference 
of the diagrams of the full theory and the corresponding diagrams in the 
effective theory. If the internal quark is heavy, logarithms of heavy masses 
will be included in the matching coefficients of the hamiltonian. If the internal 
quark is light, rrii = rric or m„, by performing the matching at fi = my/ we 
have an exact cancellation of the term |lnxi in the coefficients. In other 
words, at the matching scale all IR behavior cancels, and the coefficients are 
manifestly finite for ttij — > 0. 

2.2 RG rescaling 

Perturbative QCD corrections introduce logarithms a"(/i)ln'"(/i/mvy), with 
m < n. The RG rescales the coefficients of the effective hamiltonian to scales 
lower than the matching scale mw, and resums such logarithms. After the 
RG group rescaling, we are left with a residual dependence on the scale /i, due 
to the truncation of the perturbative series. In 5-decays, the first threshold 
is at /i = nib, and we can stop rescaling at that point; in D and K decays one 
can do a new matching and then use the RG once again to go to a still lower 



^Note, however, that the Fi term decouples from an on-shell photon with physical 
polarization. 



scale. In any case, if we are to work with perturbation theory, the final scale 
must be greater than 1 GeV or so. It is self-evident that this step cannot 
introduce logarithms of light masses. 

2.3 Matrix elements 

The final step consists in calculating the matrix elements of the operators in 
the low energy theory. In general, matrix elements include loops of virtual 
light quarks. We will study the limit where the penguin loop quark mass 
^TT-ioop goes to zero. We want to show that for the specially interesting case 
of 6 — > S7 this limit also does not introduce IR divergences. This implies 
that, in the final result, any term of the type Inrriioop will always appear 
multiphed by powers of rriioop (that is, as mf^^pln^ rriioop with a > 0). A 
number of explicit calculations corroborate this expectation; see particularly 
Ref. 0, where several matrix elements at two loops are calculated. 

Let us consider the Feynman diagrams that describe the matrix elements 
in the effective theory at arbitrary order in QCD perturbation theory. An 
arbitrary Feynman diagram with A^ lines may be written in terms of FejTiman 
parameters as 
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where Ij is the momentum of the jth line and aj its Feynman parameter {Ij is 
linear in the loop momenta kr and in the external momenta ps)- The function 
F contains overall factors that do not enter the arguments at this point. This 
integral can be viewed as a contour integral in a multidimensional complex 
aj,kr space. In order to find possible logarithms, we have to look for the 
regions of non-analyticity of the integral. The points of the contour where D 
vanishes are called "singular points" (SP's), and possible singularities of this 
integral must arise from zeros of the denominator D. In fact, only certain 
SP's, referred to as pinch SP's, give singularities that cannot be avoided by 
deforming the integration contours. Necessary conditions for a pinch SP are 
given by the so-called Landau equations |]T2| . 



With each SP is associated a reduced diagram, constructed from the com- 
plete graph by simply contracting all lines which are off-shell at the SP. The 



reduced diagrams of pinch SP's have a direct physical interpretation ^ |13 
They can be interpreted as a picture of a classical, energy- and momentum- 
conserving process occurring in space-time, with all internal particles real, on 
the mass-shell, and moving forward in time. We may turn this interpretation 
around, in order to identify SP's that are pinched. We select the reduced 
diagrams associated with an arbitrary Feynman graph that admit such a 
physical interpretation; these diagrams identify pinch SP's. 

Once we have all the reduced diagrams relevant to a particular Feynman 
graph, we know its sources of non-analyticity. At this stage, it becomes 
important to have criteria for determining which pinch SP's may actually 
introduce infrared divergences in the diagram. The presence of pinch SP's 
reveals the presence of non-analytic terms, such as logarithms of light masses. 
If these logarithms are suppressed by powers of the light masses themselves, 
however, the corresponding amplitude will not diverge in the zero-mass limit. 

In order to identify possible IR divergences we use IR power counting. 
An obvious complication for IR power counting in Minkowski space is that 
k^ = does not imply A; = 0, so that a naive dimensional counting will 
not necessarily express the real behavior of the integral in the IR limit. A 
method for dealing with this problem is to change variables, and approximate 
the integral near each pinch SP, so that every denominator is a homogeneous 



function of a set of variables that vanish there [^ |T^ . This integral will be 
referred to as the "homogeneous" integral and these variables as the "normal" 
variables; the remaining variables, called "intrinsic" , parametrize the relevant 
surface of SP's, and do not contribute directly to the singular behavior. The 
IR behavior of the homogeneous integral will be determined by dimensional 
power counting involving only the normal variables. In the following section, 
we apply the power counting procedure just sketched to b —>■ 37, and verify 
that the amplitude is free of unsuppressed logarithms of mioop for mioop — ^ 0. 



3 The decay 6 ^ 57 

The decay b ^ sj has been extensively studied in the framework of effective 



field theory [jT5[. Beyond leading order, the matrix elements of the operators 
in the effective theory include light quark loops in general. Let us consider 
one of these diagrams: precisely the penguin diagram in Fig. 2, without QCD 
corrections. There is a light quark (-u-quark or c-quark) running in the loop. 



We consider the zero mass limit for this quark. 

The pinch SP corresponding to the diagram in Fig. 2 is associated with 
a reduced diagram that coincides with the original one. At the pinch SP, 
the reduced diagram can be interpreted as a process occurring in space-time, 
with all internal particles real, on the mass-shell, and moving forward in time. 
Then it is easy to see that the two light quarks and the photon belong to the 
same jet, where a jet is defined as a connected set of massless hues, which 
are on shell with finite energy, and have momenta proportional to a single 
lightlike momentum (p^ in this case) . Therefore, the reduced diagram in Fig. 
2 can be viewed as a massive b quark decaying into two jets, one consisting 
of the light quark loop and the photon, the other consisting of the s-quark 
line only. We use the term "hard" for any vertex of a reduced diagram 
where lines from two or more jets are attached. We will also refer to on-shell 
massless lines with zero 4-momentum as soft lines; a "soft subdiagram" is 
one consisting of only soft lines. 

Let us first analyze arbitrary reduced diagrams that contribute to the 
four-quark operator with two jets, J^ and Jg, a single hard part and a soft 
subdiagram. Fig. 3 shows a typical diagrammatic structure for J^; Fig. 
4 shows the general form of these diagrams. Afterwards, we shall treat the 
remaining, relevant reduced diagrams. The diagrams of Fig. 4 admit a rather 
simple IR power counting, analogous to the treatment of form factors for 



quark- ant iquark production in e"'"e~ annihilation ||^, |TJ]. An appropriate 
choice of normal variables for all these processes is the four components 
fc^, /i = ... 3, of loops ks that pass through the soft subdiagram, S, and 
the squares k"^ and {kl _|_)^ for each loop kj that is internal to a jet. For the 
latter, the transverse momentum is defined relative to the jet's direction. 

The superficial IR degree of divergence related to any reduced diagram 
of this sort is 

D - ' 3 

={7,s} 



J2 (2L,-N, + bi + lf, + u), (3) 



where Li and Ni are the numbers of loops and lines in Jj, while hi and fi are 
the numbers of soft gluons and soft photons attached to Jj at the pinch SP. 
The factor t,- comes from the numerator momenta in J,- and is bounded from 



below0, 



U>max\-[ul-Vi],0\, (4) 



where u\ and Vi are, respectively, the number of three-point vertices in Jj 
and the number of soft vector particles attached to Jj. The Euler identity 
and counting relations between the numbers of lines and vertices of various 
orders may then be used to bound the IR degree of divergence by 



D> Y. [l{k-l) + f. + l{v,-umv.-ul)y (5) 



where hi is the number of lines attaching jet i to the decay vertex. Clearly, 
the lower limit on D is found by taking hi = 1 and fi = 0. 

We can easily check that for the lowest-order diagram. Fig. 2, D > 0, 
because in this case h^ = 2 for the photon jet, which immediately gives 
D = 1/2. We may trace the positive value oi D to ul = 1 in Eq. (^). This 
suppression is a direct result of the transversality of the emitted photon, 
which requires at least one power of the transverse momentum of the soft 
quark loop. Therefore, the diagram of Fig. 2 is IR convergent, when the 
photon is on shell. The all-order power counting expression, Eq. (j^) shows 
that this reasoning holds to any order for diagrams of the form of Fig. 4, 
since h^ = 2 for all of them when the hard vertex is a four-quark operator. 

We can readily extend this reasoning to the remaining pinch SP's that 
are relevant to 6 ^ 57. Further corrections to the four-quark operators are 
shown in Fig. 5. The only singular points involving an on-shell light quark 
penguin loop that we have not yet considered are those in which soft gluons 
attach to the b quark (Fig. 5a) and those in which the penguin loop itself is 
soft (Fig. 5b). 

First, consider additional gluons attached to the b quark. The propa- 
gator of a quark radiating a soft gluon behaves as l/2p ■ k (where k is the 
momentum of the gluon) and so contributes —1 to the power counting. For 
power counting purposes, the b quark then acts like a third jet in Eq. 0, with 
no internal loops {Li = 0), no soft external fermions {fi = 0), no numerator 
suppression (tj = 0), and with the number of its lines equal to the number 

^For this argument, we ignore unphysically-polarized gluon lines that attach the jet 
to the hard scattering in covariant gauges; they do not affect the overall power counting 
discussed here |g, |lj . 



of soft gluons attached to the b hne {Ni = hi). It is easy to verify that pinch 
SP's of this sort leave Eq. (^ unchanged. 

Second, consider the class of pinch SP's for the four-quark operator in 
which the photon is radiated by the s (or h) quark, and the light quark 
penguin loop appears as part of the soft subdiagram (Fig. 5b). This diagram 
is highly suppressed in the IR compared to those of Fig. 5a, because a fermion 
propagator with momentum k^ diverges only linearly at k^ = 0, in contrast 
to a boson propagator, which diverges quadratically. Thus, penguin loops, 
whether connected directly to the photon or not, are finite in the zero-mass 
limit. 

So far, we have restricted ourselves to light quarks in penguin loops only, 
connected directly to the operators of the effective theory. Mixing of oper- 
ators in the effective theory, however, gives rise to diagrams in which there 
is no penguin loop. We will now show that in the two-body amplitude there 
are no additional logarithms of light quark masses associated with loops in 
the photon jet, that is, collinear to the outgoing photon. 

Consider Fig. 6a, in which the gluon emerges from the hard subdiagram 
(effectively, the operator Og ~ nib sicr^'^Ta Br C^^ in the standard classifi- 
cation). These diagrams contain yet another set of pinch SP's, as shown, 
in which this gluon changes into a photon due to rescattering of a loop of 
virtual light quarks with soft gluons. We may once again use the power 
counting of Eq. (||), but this time h^ = 1, and on a diagram-by-diagram 
basis, the amplitude produces logarithms of the light quark mass. Note that 
charge conjugation invariance requires at least two soft gluons attached to 
the quark loop (two C-odd gluons cannot produce a C-odd photon), so that 
this effect appears first at three loops in the perturbative matrix element. 
Nevertheless, the contribution from any SP of this sort cancels in a gauge 
invariant set of diagrams, because the photon does not carry color. This 
result follows from the factorization of soft gluons from jets, an important 
ingredient in factorization proofs for inclusive cross sections |]TB|, |T7|. 

Leaving technicalities aside, soft gluons can couple only to the total color 
charge of a jet. In fact, the contribution of the SP of Fig. 6a may also be 
pictured as in Fig. 6b, which shows that the total effect of the soft gluons is 
to insert a nonabelian phase factor in the color product between the decay 
vertex and the gluon line. The double line represents the nonabelian phase; 



the relevant Feynman rules are described, for example, in [|T^ . Here, however. 



it is only the topology of the figure that is important. The remaining jet in 



Fig. 6b consists of a gluon and a photon, connected by the hght quark loop, 
whose color trace vanishes identically. Let us note that this factorization does 
not require a sum over final states; indeed, the cancellation of soft gluons in 



inclusive processes in Refs. [|§, [T^ requires the factorization as a first step, for 
each final state individually. Thus, logarithms of mioop associated with these 
processes are absent, since, as we have just seen, a gluon cannot fragment by 
interacting only with soft quanta. Note that this reasoning does not apply to 
final states describing the coUinear splitting of the gluon into a photon plus 
other gluons. 

In summary, we have shown, to all orders in ag and to the lowest order 
in aew, that penguin-like diagrams relevant for the process b — > s((i)7 can be 
safely calculated by taking the massless limit for the w-quark circulating in 
the penguin loop. This implies, for instance, that no large "enhancements" 
(i. e. IR unsafe contributions) in the amplitudes involving the CKM matrix 
elements VubV*^ are expected from this source at any fixed order in perturba- 
tion theory. We have also observed that in the perturbative expansion of the 
two-body final state, no logarithms of light-quark masses arise from loops 
coUinear to the outoing photon. 
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Figure Captions 

Fig. 1 One of the penguin diagrams of the full theory, as computed, for ex- 



ample, in Ref. [|TT |. 



Fig. 2 Penguin diagram of the effective (low energy) theory. 

Fig. 3 Typical reduced diagram for the QCD corrections to the penguin-like 
diagrams. 

Fig. 4 General reduced diagram for soft gluon corrections to the b -^ 2 jets 
transition. 

Fig. 5a Same as Fig. 4, but with gluonic insertions in the 6-quark line. 

Fig. 5b Reduced diagram, analogous to those of Figs. 4, 5a, but with the 7 
arising from an external (6 or s) quark line, instead of the light quark 
loop. 

Fig. 6a Reduced diagram with a gluon arising from the hard vertex. The sub- 
diagram hard gluon -^ soft gluons -|- 7 vanishes by color conservation 
(see text). 

Fig. 6b Factorization of the soft gluon contribution to the diagram of Fig. 6a. 
The double line represents the nonabelian phase referred to in the text. 
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Fig.2 





Fig.3 



Fig.4 



The hard vertex is represented by the black disk. 

The dashed lines are gluons. 

The blobs Js and Jy indicate the s and y jets respectively. 





Fig.Sa 



Fig.Sb 





Fig. 6a 



Fig.6b 



The hard vertex is represented by the black disk. 

The dashed lines are gluons. 

The blobs Js and Jy indicate the s and y jets respectively. 

The double line of Fig. 6b represents a non abelian phase. 



